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SUMMAR! 


A  method  is  presented  for  determining  the  control  gearing  and  time 
lag  necessary  for  a  specified  damping  of  the  motions  of  an  aircraft 
equipped  with  an  autopilot  having  constant-time-lag  characteristics. 

The  method  is  applied  to  a  typical  present-day  airplane  equipped  with  an 
autopilot  which  applies  rudder  control  proportional  to  the  yawing  angular 
acceleration.  The  types  of  motion  predicted  for  this  airplane -autopilot 
system  by  this  method  are  in  very  good  agreement  with  the  airplane 
motions  calculated  by  a  step-by-step  procedure. 

For  some  values  of  control  gearing  there  may  exist  more  than  one 
range  of  time  lag  for  which  the  motions  of  the  assumed  aircraft- 
autopilot  system  will  have  more  than  a  specified  amount  of  damping. 


INTRODUCTION 


Many  present-day  aircraft  designed  to  fly  at  transonic  and  super¬ 
sonic  speeds  have  exhibited  poor  lateral  stability  characteristics.  As 
a  result,  much  interest  has  been  shown  in  automatic  stabilization  systems 
as  a  means  of  obtaining  satisfactory  stability  for  high-speed  flight. 

In  analyzing  the  effect  of  a  particular  autopilot,  the  usual  practice 
has  been  to  determine  whether  the  aircraft -autopilot  system  is  stable 
by  employing  methods  such  as  those  presented  in  reference  1  (frequency- 
response  analysis)  and  reference  2  (Nyquist  plots) .  For  autopilots 
characterized  by  constant  time  lag  (linear  variation  of  control  gearing 
and  phase  lag  with  frequency)  the  same  information  can  be  obtained  by 
carrying  out  an  analysis  such  as  that  of  reference  3.  This  type  of 
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information  is  important  in  the  analysis  of  an  autopilot -aircraft 
system  hut  affords  little  q^uantitative  indication  of  the  degree 
of  stability  vhich  can  be  obtained  by  use  of  a  particular  autopilot. 

Since  one  of  the  purposes  of  equipping  an  aircraft  with  an  autopilot  is 
to  improve  its  stability^  the  determination  of  the  type  of  autopilot 
frequency  response  that  will  result  in  satisfactory  stability  appears  to 
be  of  greater  importance.  The  degree  of  stability  which  can  be  obtained 
by  equipping  an  aircraft  with  a  linear  autopilot  can  be  calculated  by 
the  methods  of  reference  4^  but^  if  an  over-all  picture  is  desired  of  the 
variation  of  stability  for  different  combinations  of  control  gearing  and 
constant  time  lag^  the  calculations  become  too  laborious  for  practical 
application.  A  method  of  obtaining  this  information  was  discussed  very 
briefly  in  the  appendix  of  reference  4,  but  has  been  found  inadequate 
for  a  comprehensive  analysis  of  autopilot  systems  which  exhibit  constant- 
control-gearing  and  constant-time-lag  characteristics. 

The  purpose  of  the  present  paper,  therefore,  is  to  extend  the 
concept  discussed  in  the  appendix  of  reference  4  and  thereby  present  a 
rigorous  method  of  obtaining  the  combinations  of  constant  control  gearing 
and  constant  time  lag  necessary  to  provide  specified  amounts  of  damping, 
up  to  the  maximum  obtainable,  to  the  motions  of  a  linear  oscillating 
system,  with  application  to  the  automatic  stabilization  of  aircraft. 

The  method  presented  is  applicable,  in  a  strict  sense,  only  to  control 
systems  which  exhibit  linear  frequency-response  characteristics,  but, 
under  certain  conditions  discussed  in  the  paper,  the  analysis  is  valid 
for  systems  having  a  frequency  response  similar  to  that  of  a  constant - 
time-lag  system  over  a  limited  range  of  frequencies. 

Although  in  this  paper  the  method  is  applied  to  the  problem  of 
automatic  stabilization  of  aircraft,  it  is,  in  general,  applicable  to 
any  linear  oscillating  system  with  constant  time  lag. 


SYMBOLS  AND  COEFFICIENTS 


0  angle  of  roll,  radians 

^  angle  of  yaw,  radians 

p  angle  of  sideslip,  radians  (v/V) 

r  yawing  angular  velocity,  radians  per  second  (d>i^/dt) 

‘•^r  yawing  angular  acceleration,  radians  per  second  per  second 

(d^^z/dt^) 
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rolling  angular  velocity^  radians  per  second  (d0/dt) 
sideslip  velocity  along  lateral  axis^  feet  per  second 
airspeed^  feet  per  second 
mass  density  of  air^  slugs  per  cutic  foot 

dynamic  pressure^  pounds  per  square  foot 

ving  span,  feet 

wing  area,  square  feet 

weight  of  airplane,  pounds 

mass  of  airplane,  slugs  (w/g);  any  integer 

acceleration  due  to  gravity,  feet  per  second  per  second 

relative-density  factor  (m/pSh) 

inclination  of  principal  longitudinal  axis  of  airplane  with 
respect  to  flight  path,  positive  when  principal  axis  is 
above  flight  path  at  nose,  degrees 

angle  of  flight  path  to  horizontal  axis,  positive  in  a  climb, 
degrees 

radius  of  gyration  in  roll  about  principal  longitudinal  axis, 
feet 

radius  of  gyration  in  yaw  about  principal  vertical  axis,  feet 

nondimensional  radius  of  gyration  in  roll  about  principal 
longitudinal  axis  ^kx^jb^ 

nondimensional  radius  of  gyration  in  yaw  about  principal 
vertical  axis 

nondimensional  radius  of  gyration  in  roll  about  longitudinal 


stability  axis  ^\/Kxq2cos2t]  +  K^Q^sin^ij^ 


nondimensional  radius  of  gyration  in  yaw  about  vertical 
stability  axis  ^^^K^q^cos^ti  +  K^^^sin^T]^ 


k 
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nondimensional  pro duct -of- inertia  parameter 
((kZo^  -  Kxo^jsin  ti  cos 


trim  lift  coefficient 


'¥  cos  7 


rolling -moment  coefficient 


yawing -moment  coefficient 


lateral-force  coefficient 


^Rollins  moment 


Yawing  moment 


Lateral  force 
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CZ.  = 


2V 


*^116 


B6r 


,  _  ^Cn 

^5a  -  ^ 


C 


25r 


SCz 

SSa 

^Cz 

BSp 

SCy 

BCy 

BSp 


t 

Sb 

I>b 

P 

^1/2 

a 

O) 

CDs 


time^  seconds 

nondimensional  time  parameter  based  on  span  (Vt/b) 
differential  operator 

period  of  oscillation^  seconds 

time  for  amplitude  of  oscillation  to  damp  to  one-half  its 
original  value ^  seconds 

deflection  of  control  surfaces^  radians 

real  part  of  complex  root  of  characteristic  stability  equation 
angular  frequency^  radians  per  second 
nondimensional  angular  frequency 

complex  root  of  characteristic  stability  equation  (a  +  ia)s) 


X 
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time  lag  between  signal  for  control  and  its  actual  motion, 
seconds 


nondimensional  time 


lag  (^t) 


amplitude  of  control-surface  oscillation  produced  by  autopilot 
in  response  to  oscillation  of  airplane  acceleration 

(h\  ) 

yl  \|f  lautopiloty 


,  =  =  vi  15. 

^  "  b2 

e  phase  angle,  radians 

EQUATIONS  OF  MOTION 

The  nondimensional  linearized  equations  of  motion,  referred  to  the 
stability  axes,  are; 

Yawing 

%Z  +  l^np  ^ 

Rolling 

+  %Z  Db^t)  =  Cjpp  +  Db0  +  (lb) 


Sideslipping 

2^b(BbP  +  I>b*)  =  +  i^fp  ^  +  (Cl  ?)♦ 

The  derivatives  Cy^^,  and  Cy^^  have  been  neglected 

in  equations  (l). 

If  a  constant-time-lag  autopilot  which  applies  rudder  control  pro¬ 
portional  to  the  nth  derivative  of  the  yawing  displacement  is  installed 
in  the  airplane,  the  equation  for  as  a  function  of  sg  becomes 


NACA  TN  23OT 


7 


/  \ 

-  Ts)  (2 

where  the  term  I>b’^\lf(st  -  Tg)  signifies  that  the  control  deflection  at 

a  given  instant  is  proportional  the  the  nth  derivative  of  t  which 
existed  at  a  fixed  time  Tg  previous  to  the  given  instant.  The  term 

o6j. 

control -gearing  ratio  of  the  autopilot.  Similarly,  if  a 

constant-time-lag  autopilot  which  applies  aileron  control  proportional 
to  the  nth  derivative  of  the  rolling  displacement  is  Installed  in  the 
airplane,  the  equation  for  69,  as  a  function  of  s^,  becomes 


^a(®'b)  “ 


ai}b°0 


Db^jZi(sb  -  Tg) 


For  purposes  of  illustration^  assume  that  the  type  of  autopilot 
described  by  equation  (2)  is  to  be  investigated.  Therefore^  the  value 


given  by  equation  (2)  is  substituted  for 


in  equations  (l). 


The  aileron  deflection 
tuted  for  0^ 


)a  is  assumed  zero.  When  0Qe 

/I  Q  b  -f*  Q  -1 


is  substi' 


for  and  PqC 


in  the  resultant 


equations  written  in  determinant  form^  X  must  be  a  root  of  the  charac¬ 
teristic  stability  eq-uation 

Ax5  +  +  C\3  +  dx2  +  ex  +  ksVe'^^HA'x3  +  B' +  C  X  +  D' )  =  0  (4) 


where  A,  B,  E,  A*,  B* ^  C’,  and  D*  are  functions  of  the  mass  and 
aerodynamic  parameters  of  the  airplane.  The  expressions  for  A^  B^  C, 
D,  and  E  are  given  in  reference  5^  and,  for  this  particular  case, 

A*  = 


=[-¥h%  ^ 


=  CLCzpCng 


and  k. 
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The  time  to  damp  and  the  period  of  the  lateral  oscillation  in  seconds 
are  given  hy  the  expressions 

T  /  -  b  1 

^1/2  -  a  V 

(5) 

P  =  2it  t 

cSs  V 

where  a  and  csg  are  the  real  and  imaginary  parts,  respectively,  of  a 
complex  root  of  equation  (4). 


AMLYSIS 


The  following  assumptions  are  made  in  the  derivation  of  the  method 
which  may  he  used  to  determine  the  combinations  of  kg  and  Tg  neces¬ 
sary  to  provide  a  specified  amount  of  damping  of  the  lateral  oscillations 
of  an  automatically  stabilized  airplane: 

(1)  The  characteristic  stability  equation  of  an  airplane  equipped 
with  an  autopilot  which  applies  rudder  control  proportional  to  the  nth_ 
derivative  of  the  yawing  displacement  which  existed  at  time  s-^  -  Tg  is 
given  by  equation  (4);  that  is,  the  autopilot  is  assumed  to  have  linear 
characteristics. 

(2)  This  stability  equation  has  as  one  of  its  roots  X,  =  a  +  iuis. 
(This  root  X  represents  an  oscillatory  mode  of  motion,  the  time  to 
damp  to  one-half  amplitude  and  the  period  of  which  are  given  by  equa¬ 
tions  (5).) 

The  stability  equation  of  the  airplane-autopilot  system  (equa¬ 
tion  (4))  may  be  rewritten  in  the  following  form: 

-Te,X  AX5  +  BX^  +  CX3  +  DX^  +  EX 

+  C').+  D') 

If  X  =  a  +  itus  is  substituted  in  each  side  of  equation  (6),  the 
condition  that  X  be  a  root  of  the  characteristic  stability  equation  is 
that  the  complex  number  A]_  +  iB]_  obtained  from  the  left-hand  side  must 

be  equal  to  the  one  obtained  from  the  right-hand  side  A2  +  iB2.  The 
quantities  Aq  +  iBq  and  A2  +  iB2  may  be  represented  by  the  expres¬ 
sions  Rqe^®!  and  respectively.  Therefore,  this  requirement 

is  equivalent  to  saying  that  Rq  =  R2  and  0q  =  ©2  ^  ® 
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root  of  equation  ( U)  .  If  X.  =  a  +  icOg  is  substituted  in  the  left-hand 
side  of  equation  (6),  the  following  expression  results: 


cos  TgCUg  -  1  sin  TgCSg 


Since  the  complex  number  A]_  +  iB]_  is  unchanged  vhen  its  phase  angle  is 
increased  by  integral  multiples  of  2^^  the  following  identities  exist: 


■1  - 


1  s  2inii  -  (Tg)  CDg 


where  m  -  0^  1^  2,  ...  .  If  X,  =  a  +  itSg  is  substituted  in  the 

right-hand  side  of  equation  (6),  the  values  of  Rg  and  9^  can  be 

calculated  nimerically.  For  each  value  of  m,  the  value  of  (t„) 

®  m 

necessary  to  make  9-^  equal  to  Sg  calculated  and  since  a  is 

a  known  fixed  quantity,  the  values  of  (k  )  necessary  to  make  R, 

°  m  1 

equal  to  R2  also  be  determined.  For  the  combinations  (ts^2_ 

(^s)-|^  (t^s)p  (^s)p^  •  •  •  ('^s^  (^s)  }  therefore, 

X  =  a  +  icDg  is  a  root  of  the  characteristic  stability  equation  (equa¬ 
tion  (4)).  The  parameter  m  in  the  preceding  equations  is  necessary 
if  the  analysis  is  to  hold  for  values  of  a  which  represent  dampings 
less  than  the  natural  airplane  damping.  More  specifically,  the  ranges 
of  Tg  and  kg  which  result  in  greater  stability  than  the  natural 

stability  of  the  airplane  alone  cannot  be  obtained  correctly  by 
utilizing  the  methods  as  presented  in  reference  4.  For  example,  the 
curves  presented  in  figure  10  of  reference  4  fail  to  provide  the 
important  information  that,  for  small  values  of  the  autopilot- 

airplane  system  discussed  in  that  reference  is  less  stable  than  the 
airplane  alone.  In  addition,  the  parameter  m,  is  necessary  if  the 

analysis  is  to  hold  for  values  of  Tg>£E-  These  values  of  Tg  must 

be  included  in  a  theoretical  analysis  of  the  problem  since,  as  T3  is 

increased  by  multiples  of  —}  the  phase  relation  again  exists  which 

s 

is  necessary  for  the  system  to  have  an  oscillatory  mode  the  damping  of 
which  is  indicated  by  a.  The  analysis  of  reference  4  considered  only 
the  case  of  m  =  1  and  was  therefore  restricted  to  values  of 

S  CUg 

or,  in  effect,  to  values  less  than  the  natural  period  of  the  oscillation. 
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For  a  given  value  of  a  corresponding  to  some  desired  value  of 
T1/2  of  the  lateral  oscillation,  the  analysis  may  he  made  throughout 

the  range  of  0)3  and  the  required  combinations  of  (ks)ni  ('^s)m 

determined  for  each  csg.  For  each  m,  therefore,  a  curve  representing 
this  value  of  a  may  then  be  plotted  as  a  function  of  (ks)jQ^  and  (rsjjn- 

This  family  of  m  curves  'will  divide  the  ks,Ts  plane  into  regions  for 
which  the  lateral  oscillations  will  have  more  or  less  damping  than  that 
indicated  by  the  value  of  a  being  investigated. 

General  expressions  for  (ts)j^  and  (ks)jj^  for  each  mg,  in  terms 
of  R2  follows: 


(Ts)m 


(ts)o 


+ 


2jtm 

^s 


(^s)tn 


where 


(■^s)o 


(ks)o  =  1^2' 


(^1 


0  =  R. 


Thus,  for  any  value  of  cjds,  it  is  necessary  only  to  obtain  (ts)o  and 
(kc.)n  since,  by  use  of  the  preceding  generating  functions,  the  required 
combination  of  Xg  and  k^  are  readily  obtained  for  any  integral  value 
of  m. 


ILLUSTRATIVE  EXAMPLES  AND  DISCUSSION 


Since  the  purpose  of  installing  an  autopilot  is  to  improve  the 
lateral  stability  of  the  airplane,  one  logical  value  of  a  to  investi¬ 
gate  would  be  that  value  which  corresponds  to  the  damping  of  the  lateral 
oscillation  of  the  airplane  without  the  autopilot  Installed,  hereinafter 
referred  to  as  airplane  damping.  Thus,  the  curves  corresponding  to  that 
value  of  a  would  indicate  those  combinations  of  kg  and  xs  foi"  which 
the  airplane  is  more  stable  with  a  given  autopilot  than  without  it.  This 
value  of  a  may  be  determined,  of  course,  by  solving  for  the  complex 
roots  of  the  characteristic  stability  equation  when  kg  is  set  equal  to 
zero  The  family  of  curves  obtained  for  this  value  of  a  forms  the 
boun^ry  between  two  different  types  of  curves;  that  is,  the  curves  for 
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dampings  greater  than  the  airplane  damping  are  of  a  different  sliape  than 
those  for  dampings  less  than  this  value.  This  change  in  the  damping 
curves  is  discussed  in  detail  in  a  subsequent  paragraph. 

For  purposes  of  illustrating  the  present  method^  the  airplane  is 
assTomed  to  have  freedom  only  in  yaw.  Several  damping  cirrves  were  calcu¬ 
lated  for  a  typical  present-day  airplane  equipped  with  an  autopilot 
sensitive  to  yawing  acceleration  for  both  the  three -degree -of -freedom 
and  the  one -degree -of -freedom  cases.  The  results  of  the  calculations^ 
although  not  always  in  quantitative  -agreement,  indicated  that  the  nature 
of  the  damping  curves  obtained  for  both  cases  was  essentially  the  same 
and  therefore  subsequent  discussion  is  based  on  the  results  obtained  when 
the  system  was  assumed  to  have  freedom  only  in  yaw. 


The  mass  and  aerodynamic  characteristics  of  the  typical  present-day 
airplane  are  presented  in  table  I.  The  characteristic  stability  equa¬ 
tion  for  the  three-degree-of -freedom  case  was  presented  previously 
(equation  (4)),  and  the  characteristic  stability  equation  for  the  one- 
degree  -of  -freedom  case  is  given  by 

(7) 


The  derivative  is  assumed  equal  to  -Cnj3.  The  equation  used  in  the 

calculations,  obtained  by  rearranging  equation  {!) ,  is 


kge 


-  iCn^X  - 


(8) 


As  was  pointed  out  in  a  previous  section,  the  damping  curve  corre¬ 
sponding  to  the  damping  of  the  airplane  oscillatory  motion  without  the 
autopilot  installed  forms  the  boundary  between  two  characteristically 
different  types  of  curves.  In  order  to  illustrate  the  types  of  curves 
which  will  be  encountered,  calculations  were  made  for  the  airplane 
damping  and  for  dampings  less  than  and  greater  than  this  damping.  The 
results  are  shown  in  figure  1.  For  the  airplane  described  in  table  I, 
the  airplane  damping  is  Tp/p  ”  2.02  seconds.  The  values  of  Tg  and  kg 
obtained  from  equation  (8)  have  been  converted  from  the  nondimens ional 
time  parameter  Sp,  to  time  in  seconds.  Thus,  the  curves  are  plotted  as 
a  function  of  the  dimensional  parameters  k  and  t  in  seconds.  Each 
point  on  these  curves  corresponds  to  a  different  frequency  cOg,  which 
varies  from  infinity  at  the  point  (T=0,k=0.0628)  to  zero  as  t  and  k 
approach  infinity.  It  is  interesting  to  note  that  all  the  damping 
curves  approach  the  point  (0,0.0628)  as  (Dg oo.  For  this  point,  the 
roots  of  the  characteristic  stability  equation  do  not  indicate  any  high- 
frequency  oscillation;  however,  for  any  finite  value  of  t,  the  roots  of 
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this  equation  will  indicate  a  high-frequency  oscillation  in  addition  to 
the  usual  lateral  oscillation.  For  a  given  value  of  k,  the  frequency 
of  this  additional  oscillation  decreases  with  increases  in  t.  Also, 
all  the  damping  c\arves  approach  a  minimum  value  of  k  for  a  value  of  cOg 
corresponding  to  a  frequency  close  to  the  natural  frequency  of  the  air¬ 
plane.  In  general,  frequencies  greater  than  and  less  than  the  airplane 
natural  frequency  are  separated  "by  the  point  of  minimum  k  on  each  of 
the  damping  curves.  For  this  airplane  the  natural  frequency  is  approxi¬ 
mately  5  radians  per  second. 

The  curves  presented  in  figure  l(a)  correspond  to  ^ 

complete  boundary  for  neutral  stability  is  made  up  of  the  curves  for 
m  =  0,  1,  2,.  .  .  “.  Unless  m  is  allowed  to  take  on  integral  values 
it  is  impossible  to  define  correctly  the  stable  region  of  the  k,T  plane. 
For  example,  consider  the  point  (k=0.035,'r=l.  6)  in  figure  l(a).  The 
point,  according  to  the  cirrve  for  m  =  1,  represents  a  stable  condition, 
but  since  it  falls  on  the  unstable  side  of  the  m  =  2  curve,  the  pre¬ 
sence  of  unstable  oscillatory  modes  in  addition  to  the  stable  oscillation 
is  indicated  and  thus  this  combination  of  k  and  t  would  result  in  an 
unstable  condition.  The  true  boundary  of  neutral  stability  is  shown  as 
the  solid  curve  of  figure  l(a).  Only  the  curves  for  m  =  1,  2,  and  3 
are  shown  The  curve  for  m  —  0  would  fall  in  the  region  corresponding 
to  negative  values  of  t  and  therefore  was  omitted.  From  this  figure 
it  is  apparent  that  for  some  values  of  k  there  exist  several  ranges 
of  T  for  which  the  system  will  be  stable.  The  width  of  these  ranges 
diminishes  as  m  increases  and  eventually  disappears  completely.  For 
this  special  case  of  T1/2  =  »  and  a  given  value  of  k,  the  ranges 
of  T  necessary  for  stability  can  be  calculated  analytically  by  the 
methods  of  reference  3^  s-ince  the  airplane  was  originally  stable.  The 
results  obtained  by  using  this  method  agree  exactly  with  those  of 
figure  l(a).  The  analysis  presented  in  this  reference  was  for  a  one- 
degree-of -freedom  system  and  in  appendix  A  it  has  been  modified  to  apply 
to  an  airplane -autopilot  system  with  three  degrees  of  freedom.  It  should 
be  emphasized  that  the  method  of  reference  3  is  applicable  only  to  the 
special  case  for  T^^/g  “  "  originally  stable  airplane.  On  the 

other  hand  the  method  presented  in  this  paper  is  equally  applicable  to 
originally  unstable  airplanes,  in  which  cases  the  T^/g  =  curves  will 

have  the  form  of  the  curves  for  dampings  greater  than  the  original  air¬ 
plane  damping. 

For  '^±/2  ^  3*5®  seconds,  the  resulting  curves  for  m  =  0,  1,  2, 

and  3  are  presented  in  figure  l(b) .  The  general  shape  of  the  curves  is 
the  same  as  for  those  of  figure  l(a),  and  the  same  general  conclusions 
apply.  Only  that  part  of  the  curve  for  m  =  0  which  is  in  the  region 
of  positive  time  lags  is  shown  in  the  figure.  The  curves  for  the  air¬ 
plane  damping,  T1/2  =  2.02  seconds,  are  shown  in  figure  l(c).  For  the 
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combinations  of  k  and  t  within  the  triangular -like  regions,  which 
are  formed  by  the  curves  for  m,  m  +  1,  and  k  =  0,  the  oscillatory 
motion  of  the  airplane  is  more  stable  with  the  autopilot  than  without 
It.  The  damping  curves  for  T1/2  =  1.4o  seconds  and  Tp/g  =  O.70  second 
are  presented  in  figiires  1(d)  and  1(e),  respectively,  for  m  =  1  2 

and  3.  These  curves  are  quite  different  from  the  ones  presented'for 
ping  less  than  the  airplane  damping.  The  combinations  of  k  and  t 
for  which  the  airplane  would  be  more  stable  than  is  indicated  by  the 
curves  for  T1/2  =  1.^  seconds  or  T1/2  =  0.70  second,  for  example, 
are  enclosed  in  island-like  regions  in  the  stable  part  of  the  k,T  plane 
A  similar  region  exists  for  each  value  of  m,  but  as  m  is  increased 
the  size  of  the  region  diminishes.  In  addition,  as  Tp/2  is  decreased 
these  looped  regions  become  increasingly  smaller  for  each  value  of  m 
As  T1/2  is  reduced  beyond  some  finite  value  these  loops  will  cease  to 
exist.  The  value  of  T1/2  ^or  which  the  m  =  1  loop  disappears 


represents  the  maximum  damping  obtainable  for  this  airplane  by  the  use 
o  a  specific  autopilot.  A  method  of  obtaining  this  maximum  damping 
for  each  value  of  m  is  presented  in  appendix  B.  It  can  be  seen  f^om 
igure  1(e)  that  no  loops  exist  for  m  >  1.  This  amount  of  damping, 
therefore,  can  only  be  obtained  for  all  the  oscillatory  modes  within 
the  m  -  1  loop.  This  conclusion  could  have  been  reached  by  carrying 
out  the  analysis  of  appendix  B.  As  a  matter  of  Interest,  for  this  parti¬ 
cular  autopilot-airplane  system,  the  m  =  1  loop  ceases  to  exist  for 
T1/2  <  0.35  second. 


The  regions  in  the  k,T  plane,  which  correspond  to  the  dampings 
shown  in  figure  1,  have  been  combined  and  are  presented  in  figure  2. 

This  figure  gives  an  over-all  pictirre  of  the  effects  of  the  specific 
autopilot  bn  the  lateral  stability  of  the  airplane.  In  order  to  verify 
the  results  indicated  by  this  figure,  airplane  motions  were  calculated 
for  the  combinations  of  k  and  t  identified  as  points  A,  B,  C,  and  D. 
Step-by-step  calculations  were  made  by  ass-uming  freedom  only  in  yav,  and 
the  airplane  motions  in  yaw  subsequent  to  a  5^  displacement  in  yaw  for 
these  combinations  of  k  and  x  are  presented  in  figure  3.  Point  A 
is  located  within  the  m  =  1  loop  corresponding  to  Tj_y2  =  seconds 

and  the  motion  of  flg\are  3(a)  checks  this  prediction.  The  motion  for 
point  B  (fig.  3(li))  is  very  lightly  damped  as  is  indicated  by  its  near¬ 
ness  to  the  zero-damping  curve.  Point  C  represents  a  combination  of  k 
and  T  for  which  the  airplane  is  less  stable  with  this  autopilot  than 
without  it.  Also,  as  seen  from  figure  2,  the  motion  for  this  combina¬ 
tion  of  k  and  x  should  indicate  two  oscillatory  modes  of  different 
frequency  with  approximately  the  same  amount  of  damping.  The  motions 
of  figure  3(c)  support  this  prediction  since  more  than  one  oscillatory 
mode  is  indicated.  The  combination  of  k  and  x  for  point  D  is  located 
within  the  m  =  2  loop,  T1/2  ^  l.ij-0  seconds.  The  motion  for  this  combi¬ 
nation  (fig.  3(<i))  is  approximately  as  heavily  damped  as  that  of 


14 


NACA  TN  2307 


figure  3(a) j  but  because  of  the  larger  value  of  time  lag  this  high 
damping  is  not  realized  as  quickly.  It  is  of  significance,  however,  that 
this  amount  of  damping  can  be  obtained  for  the  combinations  of  k  and 
T  within  the  m  =  2,  3j  •  •  •  loops  as  well  as  for  those  within  the 
m  =  1  loop.  In  actual  practice  the  phase  lags  encoimtered  in  conventional 
autopilots  which  can  be  represented  by  second-order  systems  generally  do 
not  exceed  I8OO  at  any  frequency  and  the  assumption  of  constant-time-lag 
characteristics  usually  is  valid  only  within  the  frequency  range  where 
the  autopilot  phase  lag  is  less  than  approximately  960,  The  time  lags 
common  to  the  m  =  2,  3,  •  •  •  loops  shifts  greater  than 

2jt  and  thus  are  beyond  the  range  where  the  assumption  of  constant  time 
lag  is  usually  valid.  Two  conceivable  cases  may  be  mentioned  for  which 
it  might  be  desirable  to  develop  devices  with  such  large  time  lags. 

One  is  the  case  of  a  high-natural-frequency  system  (not  necessarily  an 
aerodynamic  system) ,  for  which  it  might  be  impractical  to  obtain  the 
small  lags  necessary  for  the  stability  of  the  system  with  automatic 
stabilization  in  the  m  =  1  region.  Another  case  is  that  of  an  aero¬ 
dynamic  system  which  is  to  respond  to  automatic  control  in  addition  o 
being  automatically  stabilized,  in  which  case  it  might  be  desirable _to 
delay  the  operation  of  the  stabilization  device  in  order  to  reduce  its 
effect  on  the  response  time  of  the  automatic  control. 

It  should  be  reemphasized  that  the  preceding  analysis  was  made  for 
an  airplane  equipped  with  an  autopilot  having  constant  time  lag.  Damping 
curves  such  as  "those  presen.'ted  in  figure  2  "therefore  apply^  in  a  s'tric 
sense,  only  to  autopilots  which  have  these  assumed  characteristics.  If, 
however,  a  frequency-response  analysis  of  a  given  autopilot  indica"tes 
that  the  assumption  of  a  constant  time  lag  is  valid  for  only  a  limited 
range  of  frequencies,  the  preceding  methods  can  still  be  used  to  good 
advantage.  Damping  curves  such  as  those  of  figure  2  can  be  construe  e 
and  the  values  of  k  and  t,  which  will  result  in  a  given  amount  of 
damping  of  the  oscillatory  modes  of  motion,  can  be  determined  since  for 
each  point  in  the  k,T  plane  there  will  be  certain  complex  solutions  of 
the  airplane-autopilot  characteristic  stability  equation  (equation  w J ; • 
If  the  frequencies  of  the  oscillatory  modes  characteristic  of  a  parti¬ 
cular  point  in  the  k,T  plane  are  such  that  the  autopilot  can  be  made  to 
exhibit  constant-lag  characteristics  up  to  and  beyond  the  highest 
frequency  involved,  then  the  damping  of  the  oscillatory  motions  for  the 
airplane-autopilot  system  will  be  essentially  the  same  as  predicted  by 
the  constant-time-lag  analysis.  The  validity  of  this  statement  can  be 
seen  by  examination' of  equation  (8).  The  left -h^d  side  of  this  equation 
is  the  autopilot  frequency-response  function  6r/^  Likewise  the 

right-hand  side  of  the  equation  is  the  airplane  frequency-response  func- 
tion  bp/D^2^.  necessary  and  sufficient  condition  that  X  -  a  +  loig 

be  a  solution  of  the  characteristic  equation  (equation  (T))  is  that, 
when  a  +  icOg  is  substituted  for  X,  in  equation  (8) ,  the  complex  number 
obtained  from  the  autopilot  frequency-response  function  5r/l>b^^ 
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identical  with  that  obtained  from  the  airplane  frequency-reeponee  func¬ 
tion  frequencies  characteristic  of  any  point  in  the  k,T 

plane  are  readily  identifiable  and,  if  the  frequency-response  function 
6r/l>b2t  of  an  actual  autopilot  can  "be  closely  approximated  "by  the  ex- 

-  X 

preeslon  =  k  e  ^  up  to  and  beyond  those  frequencies,  then  the 

complex  roots  of  the  airplane -autopilot  characteristic  equation  must 
be  essentially  the  same  as  those  predicted  by  a  constant- time- lag 
analysis. 


CONCLUDING  REMARKS 


A  method  Is  presented  for  determining  the  control  gearing  and  time 
lag  necessary  for  a  specified  damping  of  an  aircraft  equipped  with  a 
constant-time-lag  autopilot.  The  method  Is  applied  to  a  typical  present- 
day  airplane  equipped  with  an  autopilot  which  applies  rudder  control 
proportional  to  the  yawing  angular  acceleration.  The  results  calculated 
for  an  airplane-autopilot  system  by  the  method  described  are  In  excellent 
agreement  with  the  airplane  motions  calculated  by  a  step-by-step  procedure. 

The  Investigation  shows  that,  for  some  values  of  control  gearing, 
there  exists  more  than  one  range  of  time  lag  for  which  the  assumed 
aircraft -autopilot  system  will  have  a  specified  amount  of  damping.  The 
width  of  these  ranges  diminishes  with  Increased  time  lags,  and  beyond 
some  value  of  lag,  these  ranges  cease  bo  exist. 

It  should  be  emphasized  that  the  analysis  presented  In  this  paper 
Is  applicable,  in  a  strict  sense,  only  to  systems  which  exhibit  constant- 
time-lag  characteristics;  however,  under  certain  conditions  discussed 
in  the  paper,  a  constant -time -lag  type  of  analysis  is  valid  for  systems 
which  have  constant  time  lags  for  only  a  limited  range  of  frequencies. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va,  ,  September  26,  I95O 
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APPENDIX  A 


EXTENSION  OF  ANSOFF'S  STABILITY  CRITERION  TO  THE  LATERAL 
MOTIONS  OF  AN  AIRPLANE  EQUIPPED  WITH  YAW  AUTOPILOTS 


Reference  3  analyzes  the  effects  of  constant  time  lag  in  the  feed- 
hack  of  a  one-degree-of -freedom  oscillatory  system  on  the  stability  of 
the  oscillation.  The  analysis  shows  that^  in  general,  for  a  given 
gearing,  ranges  of  time  lag  will  exist  for  which  the  system  remains 
stable.  These  ranges  are  separated  from  each  other  by  unstable  ranges. 
The  expression  for  these  stable  time-lag  ranges  is  obtained  in  closed 
form  in  terms  of  the  other  parameters  of  the  system.  It  should  be 
emphasized  that  the  analysis  is  valid  only  for  a  system  originally 
stable  without  feedback. 

This  appendix  extends  the  method  of  reference  3  to  the  three  degrees 
of  freedom  of  lateral  motion  of  an  airplane  with  a  yaw  autopilot.  The 
argument  and  analysis  presented  in  reference  3  are  repeated  herein  only 
to  the  extent  necessary  to  present  the  modifications  of  the  method  and 
to  emphasize  or  clarify  points. 

The  autopilot  is  assumed  to  provide  a  yawipg  moment  proportional 
to  the  yaw  displacement  or  to  the  nth  time  derivative  of  the  yaw,  with 
the  moment  lagging  the  disturbance  by  a  constant  time^  t.  resulting 

nondimenslonal  equations  of  motion  are  given  by  equations  (1)  and  (2), 
and  the  characteristic  stability  equation  of  the  system  is  equation _( A) . 
Note  that  the  first  five  terms  of  equation  (A)  are  the  characteristic 
function  of  the  airplane  without  the  autopilot  installed.  If,  therefore 
equation  (A)  is  divided  through  by  this  characteristic  function  of  the 
airplane  alone,  the  characteristic  equation  of  the  total  system  is 
obtained  in  the  form 


F(X)  =  1  + 


(A'x3  +  B*x2  +  C'X  +  D’)  X^kse 


AX'’  +  BX 


^  +  cx3  +  Dx2  +  EX 


=  1  +  W(X)  =  0 


The  theory  may  be  given  very  briefly  as  follows.  The  roots  of 
equation  (A)  are  also  the  roots  of  equation  (A1)  and  are  the  character¬ 
istic  roots  of  the  system.  In  order  for  the  system  to  be  stable,  their 
real  parts  must  all  be  negative.  In  the  complex  X-plane  a  closed  curve 
consisting  of  the  imaginary  axis  and  the  infinite  semicircle  to  its 
right  will  enclose  all  values  of  X  with  positive  real  parts.  Then  by 
a  theorem  in  the  theory  of  functions  of  a  complex  variable  (see  refer¬ 
ence  6),  it  can  be  shown  that  the  transformation  of  this  closed  curve 
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into  the  complex  F-plane  will  circle  the  origin  in  the  F-plane  once  in 
a  positive  direction  for  each  root  of  F(  X,)  =0  enclosed  hy  the  curve 
in  the  X.-plane  and  once  in  the  negative  direction  for  each  pole.  It  can 
he  seen  from  equation  (a1)  that  the  origin  of  the  F-plane  corresponds 
to  the  point  (-l,0i)  in  the  W-plane,  and,  in  fact,  each  ciirve  in  the 
F-plane  is  simply  displaced  unit  distance  to  the  left  in  the  W-plane. 

If  the  airplane  were  originally  stable,  F(  X,)  would  have  no  poles  in 
the  right  half  of  the  X.-plane;  therefore,  the  condition  that  the  whole 
system  have  no  roots  in  this  right  half  of  the  plane  (unstable  roots) 
is  simply  that  W(X)  should  not  enclose  the  (-l,0i)  point  in  the 
W-plane. 

One  condition  for  the  preceding  theorem  to  apply  is  that  there  be 
no  zeros  or  poles  of  F( X)  on  the  boundary  curve  in  the  X-plane.  It 
will  be  shown  subsequently  that  for  n  >  2  stability  is  impossible; 
thus,  the  only  autopilots  of  interest  are  those  where  n  =  0,  1,  and  2, 
that  is,  displacement,  velocity,  or  acceleration  autopilots.  Now  the 
zero  root  of  the  characteristic  equation  (if)  for  n  =  1  and  2  is 
removed  by  the  division  which  gives  equation  (A1);  therefore,  no  diffi¬ 
culty  occurs  for  these  two  cases.  For  n  =  0,  however,  equation  (Al) 
gives  a  pole  for  F( X)  at  X  =  0,  which  is  on  the  boundary  curve; 
therefore,  for  n  =  0  dividing  through  by  X  is  not  convenient.  Thus, 
for  n  =  0 

’x3  +  B'x2  +  C’X  +  DOkse"”^®^ 

- 5 - ^ - 5 - 2 -  =  0  (A2) 

AX  +  BX^  +  CX"^  +  DX  +  E 


F(  X)  =  X  +  W(  X)  =  X  + 


(a 


Finally,  it  is  shown  in  reference  3  that,  because  of  the  negative 
exponential  in  equations  (Al)  and  (A2),  the  transform  in  the  W-plane  of 
the  boundary  cfirve  of  the  right  half  of  the  X-plane  degenerates  to  simply 
the  transform  of  the  Imaginary  axis.  Thus,  this  curve  can  be  obtained 
by  placing  icog  for  X  in  equations  (Al)  and  (A2)  and  letting  cOg 
vary  from  -00  to  00,  Equations  (Al)  and  (A2)  become 

F(icDg)  =  1  +  W(ia)s)  =0  (A3) 

for  acceleration  and  velocity  autopilots,  and 

F(icDs)  =  icog  +  W(lc0s)  =  0  (Alf) 

for  displacement  autopilots. 

For  displacement  autopilots,  the  critical  point  in  the  W-plane 
corresponding  to  the  origin  in  the  F-plane  becomes  the  point  -itDg,  a 
varying  point;  that  is,  W(oog)  must  not  enclose  the  point  (0,-i(Dg). 
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Acceleration  Autopilots 

The  expression  for  W  as  a  function  of  cus  ins-y  he  obtained  for 
the  acceleration  autopilot  hy  use  of  eqmtions  (Al)  and  (A3).  For  the 
acceleration  autopilot,  letting  a  =  A'kg,  h  =  B'kg,  c  =  Ckg,  and 
d  =  D'ks  gives 

W(a)s) 


and 


(ATa) 

(A7b) 


The  angles  and  0^  i^^y  he  obtained  from  these  equations  for  any 

given  value  of  eng  since  the  algebraic  sign  of  the  numerator  of  tan  9 
is  the  algebraic  sign  of  sin  0  and  the  algebraic  sign  of  the  denomi¬ 
nator  is  that  of  cos  0.  These  signs  determine  the  proper  quadrant 
for  0. 

The  W(cos)  curve  cannot  encircle  the  point  (-l,Oi)  in  the  W-plane 
unless  |w|  >1.  Thus,  the  ranges  of  a)g  must  be  fo^md  where  |w|  >  1 
by  setting  |w|  =  1  and  solving  for  eng.  In  the  following  discussion 
only  the  part  of  the  curve  for  positive  ojg  will  be  considered  inasmuch 
as  the  curve  for  negative  uig  is  simply  the  mirror  image  of  the  curve 
for  positive  (X>s  “the  real  axis. 

An  important  result  can  be  obtained  by  considering  what  happens 
to  W(a)s)  as  Og oo.  The  phase  angle  of  W  is  seen  to  spiral  around 
the  origin  whereas 
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lim  W 

cjDc;  co' 


a 

A 


kgCng 

.-.—  —  - - -  ...  ■  ■  .  ...  . . 

■  %z) 


(a8) 


Thus^  W(a)s)  spirals  about  a  circle  of  this  radius^  and  it  can  be  seen 
that  stability  is  impossible  if  the  limit  is  greater  than  1,  since  the 
point  (-l,Oi)  would  certainly  be  enclosed  in  this  case.  Therefore^ 
stability  is  impossible  for  n  >  2  since  in  this  case  lim  |w|  =  «>, 

In  the  cases  of  interest  the  limit  of  equation  (a8)  is  less  than  1. 
The  next  step  is  to  find  the  critical  values  of  CDg  by  setting  |W  |  =  1. 
This  condition  gives  a  quartic  in  ^  substituted  for 

the  critical  roots  are  the  positive  real  roots  of  the  resulting  eqioation 

(a^  -  A^)x^  4-  (b^  -  +  2AC  -  2ac)x3  +  (c^  -  +  2BD  -  2bd  --  2Ae)x^  + 


(d2  -  +  2Ce)x  -  e2  =  0 


(A9) 


When  there  are  no  positive  real  roots  of  equation  (A9)^  the  system 

is  always  stable.  Moreover,  for  jj|  <  1,  there  must  always  be  an  even 

number  of  such  roots  (double  roots  being  counted  twice).  For  the  one- 
degree-of-freedom  case  considered  in  reference  3  this  equation  was  a 
quadratic;  consequently,  there  was  only  one  pair  of  values  ^0332^^82) 

which  |w|  =1,  and  between  these  values  |w|  >'l.  This  condition  corre¬ 
sponds  to  the  fact  that  the  frequency-response  curve  of  the  system  has 
only  one  maximum.  In  the  three -degree -of -freedom  case  there  are  four 
possible  positive  real  roots,  a  condition  which  corresponds  to  the  fact 
that  the  frequency-response  curve  may  have  two  maximums  (see  fig,  4(a)). 
Since  either  of  these  maximums,  in  case  four  critical  roots  exist,  may 
enclose  the  point  (-l,Oi),  each  pair  of  successive  roots  must  be  treated 
in  the  same  way  as  the  single  pair  was  treated  in  reference  3.  In  most 
cases  equation  (A9)  will  have  only  two  positive  real  roots.  A  necessary 
condition  that  there  be  four  positive  roots  is  that  equation  (A9)  have 
four  sign  changes  in  its  coefficients.  The  necessary  and  sufficient 
conditions  for  existence  of  four  real  roots  can  be  found  in  most  algebra 
books  (for  example,  reference  7). 

Although  equation  (A9)  will  usually  have  only  one  pair  of  real 
positive  roots,  for  the  sake  of  generality  the  possibility  of  two  pairs 
must  be  considered.  Corresponding  to  the  two  pairs  of  successive  roots 
and  and  (cug^’  and  where  CDg^’  >  cOg^'  >  cug^  >  cjQg^, 

two  pairs  of  values  (  02^  (^1*^  i^2*  )  obtained  from  equa- 

tions  (A6)  and  (AT).  Then,  there  are  two  conditions  to  determine  the 
stable  ranges  of  time  lag: 
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(2m  -  l)jr  +  02  ^  ^  ^  +  l):t  +  ^ 


(AlOa) 


O), 


S2 


CO, 


SI 


(2m-l)7r  +  0p*  (2m  +  l)jt  +  0t' 

- ^ <  T  < - : - - 


(AlOb) 


CO. 


COc 


For  the  system  to  be  stable^  both  of  these  must  be  satisfied;  that  is^ 
the  stable  time-lag  ranges  are  those  where  the  ranges  in  equation  (AlOa) 
and  equation  (AlOb)  overlap.  Here  m  =  0,  1,  2,  .  .  .  .  The  values 
of  m  for  which  the  stable  ranges  exist  are  given  by  the  condition 


m  < 


(jt  +  0]^)a)s2  +  ~  02)^3^ 

2jt^a)s^  -  "sg) 


(Alla) 


or 

«  +  jZil')u)s2'  +  («  -  02 
2«((i>si'  -  ‘^2') 

Of  course  the  smaller  of  these  two  bounding  values  is  actually  the 
largest  m  for  which  stable  ranges  exist.  This  value  can  be  used  in 
the  corresponding  eqiaation  (AlOa)  or  (AlOb)  to  find  the  maximum  value 
of  T  for  stability. 


m  < 


(Allb) 


Velocity  Autopilots 


The  procedure  for  velocity  autopilots  is  the  same  as  that  for 
acceleration  autopilots,  the  only  difference  being  that  the  magnitude 
of  W(a)s)  is  divided  by  cug  and  the  phase  angle  decreased  by  jt/2. 
The  change  in  phase  angle  appears  in  01,  which  is  now  given  by 


ccOc  —  acOc^ 
tan  ei  =  — ^ ^ 

d  -  bcDg^ 


(A12) 


Setting  |w|  equal  to  1  now  yields,  for  =  x, 

+  ^b2  -  a2  -  2Ac)x3  +  (c2  -  b2  +  2ac  +  2AE  -  2Bd)x2  + 


(d^  -  c2  +  2bd  -  2Ce)x  +  -  d^  =  0 


(AI3) 
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It  might  be  noted  that  for  this  case  lim  w|  =0.  so  that  no  condition 

COg-^Ool 

similar  to  that  of  eqmtion  (a8)  arises.  However,  in  this  case,  |w| 
does  not  vanish  as  cDg  ->0: 


lim„  W  = 


CnrCZp  -  Cz^Cnp  +  -  CnpCZp^tan  y 


(Al4) 


Since  |E |  is  generally  very  small,  this  quantity  will  generally  be 
greater  than  1.  In  this  case,  equation  (A13)  will  have  one  or  three 
critical  roots  (usually  one),  as  may  be  seen  from  figures  4(b)  and  4(c). 

In  the  case  |^|  <1  an  even  number  of  critical  roots  will  exist,  and 

this  case  is  treated  exactly  as  in  the  case  of  the  acceleration  autopilot 
by  using  equations  (AlO)  and  the  associated  discussion,  with  equa¬ 
tion  (AI2)  replacing  equation  (A7a).  However,  in  the  more  general  case 

of  l^j  >  1,  a  slightly  different  treatment  is  needed. 


For  the  usual  case  of  an  odd  number  of  critical  values  of  cug,  the 
lowest  value  is  treated  separately.  Call  this  value  Then,  |w|  >1 

when  CDg  <  uisq,  as  can  be  seen  in  figure  4(b) .  Since  the  phase  angle 
(-<nsT  +  0)  is  a  monotonically  decreasing  function  of  £%  and  vanishes 
in  this  case  for  cug  =  0,  the  first  condition  for  stability  becomes 

“^so’’’  00  ^  where  0o  is  obtained  when  cngQ  is  used  in  equa¬ 

tions  (AI2),  (ATb),  and  (a6).  Thus,  the  stability  condition  for  the 
first  maximum  of  |W |  is  the  stringent  one  that 


T  < 


+  00 

“so 


(AI5) 


In  the  case  of  one  root,  this  is  the  only  condition.  For  three  roots  an 
additional  condition  must  be  satisfied  since  the  second  maximum  may 
enclose  the  point  (-l,0i).  In  this  case  the  remaining  two  roots,  with 
‘^Sq  <  “S2  <  “s^^  paired  just  as  in  the  acceleration  autopilot,  except 

that  equation  (A12)  of  course  replaces  (A7a),  and  eqviatlon  (AlOa)  applies. 
Note  that  both  conditions  (equations  (AlOa)  and  (AI5))  must  be  satisfied 
for  stability,  as  described  in  the  discussion  of  the  acceleration  autopilot. 


Displacement  Autopilots 

In  comparing  equation  (A2)  with  equation  (Al)  for  n  =  1,  the  W(ais) 
in  equation  (a4)  for  displacement  autopilots  can  be  seen  to  be  exactly  the 
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same  as  the  W(a)s)  in  equation  (A3)  for  velocity  autopilots.  The 
difference  hetween  the  two  cases  arises  from  the  change  in  the  critical 
point  of  the  W-plane,  as  pointed  out  in  the  paragraph  following  equa¬ 
tion  (a4)  .  Only  positive  values  of  0)3  need  he  considered,  and  for 
these  values  the  critical  point  is  (0,-in)),  Since  this  point  cannot  be 
enclosed  except  when  |w|  >  csg,  the  critical  values  of  are  given 

l^y  -t;]2e  positive  real  roots  of  ~  and  the  critical  phase  angles 

are  now  -^2m  +  rather  than  -(2m  -  1)tc. 

The  equation  determining  the  critical  values  for  x  =  0)3  is 

A^x^  +  (b2  -  2Ac)x^  +  (c2  -  a2  +  2AE  -  2Bd)x3  + 

(d2  _  +  2ac  -  2Ce)x2  +  (e^  -  +  2bd)x  -  d^  =  0  (Al6) 


This  equation  gives  the  intersections  of  the  function  |w|  with  the 
straight  line  with  unit  slope,  as  shown  in  figures  4(d)  and  4(e).  From 
these  figures  the  equation  is  seen  always  to  have  an  odd  number  of  real 
positive  roots  (double  roots  being  counted  twice);  and^  in  particular^ 
figure  4(e)  shows  how  the  maximm  of  five  critical  values  might  occur. 


The  method  of  procedure  is  the  same  as  that  given  for  the  velocity 
autopilots  with  an  odd  niomber  of  roots.  Equations  (A12),  (A7b) ,  and 
(a6)  are  used  to  get  0(cns)*  The  first  critical  frequency  ojsq  gives  as 

the  stability  condition  in  this  case  that  -cssqt  +  (ZIq  >  “f*  Thus,  the 

first  stability  condition  for  this  case  is  apparently  the  very  stringent 


one  that 

1+00 

T  <  ^ - 

“so 


(AIT) 


For  three  critical  values,  the  next  two  are  paired  as  before  and 
give  the  additional  stability  condition  for  cUsq  <  <ns2  “s^* 


which  replaces  equation  (AlOa).  Eqmtions  (AIT)  and  (A18)  must  both  be 
satisfied  for  stability. 


Finally,  for  five  critical  values,  an  almost  impossible  case,  the 
final  pair  of  values  gives  a  third  necessary  condition  when 

used  in  equation  (A18) ,  as  described  in  the  discussion  of  equation  (AlOb). 
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APPENDIX  B 

DETERMINATION  OF  MAXIMUM  DAMPING  OBTAINABLE  BY  USE 
OF  A  SPECIFIC  AUTOPILOT 


The  conditions  which  must  he  satisfied  if  the  damping  curves  are  to 
have  loops  are  that  ('’'s)^  have  Doth  a  maximum  and  a  minimum  value  for 
some  values  of  say  (cos)i  and  ((ns)2>  and  that  (ks)jjj  be  a  mini¬ 

mum  for  some  value  of  0)3  which  is  between  (cDs)i  and  (c0s)2-  (See 

fig.  1(d).)  The  loops  will  break  down  when  — V  and  ^ vanish 

dais  aa)s 

for  the  same  value  of  csg.  The  value  of  a  for  which  these  conditions 
are  satisfied  is  the  maximum  damping  of  the  system  obtainable  by  use  of 
a  specific  autopilot.  A  curve  representing  a  damping  greater  than  this 
maximim  in  the  kg^r  g  plane  will  have  no  practical  significance,  since 
for  every  point  on  this  c\irve  some  oscillatory  mode  with  less  damping 
than  the  maximum  will  always  exist. 

Now,  (ks)jji  and  (ts)iii  are  defined  by  the  expressions 


(Bl) 


where  R2,  (ts)]ji,  and  Og  are  functions  of  cug  and  were  discussed  in 

the  section  entitled  "Analysis.”  Differentiation  of  equations  (Bl) 
results  in  the  expressions 


d-(j^s)m 

dn)s 


[(Ts)ma] 


(Ts)ma  dR2 


<i(^s)m 

doig 


dep 

-cog  —  .  2m  +  62 

- 


Setting  these  expressions  equal  to  zero  and  considering  m  =  1  gives 


2k 
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=  aRp  — ^  ^  =  0 

ctos  <tos 


dT  d.0p 

^  =  2jt  -  ^52  +  “s  :r-^  =  0 

(tog  ®  diDg 


Since  =  0  is  a  necessary  condition,  ^  must  also  be  equal  to 

d£Dg 

zero  if  is  to  vanish.  The  expressions  for  R2  and  @2 

d£i)g 

functions  of  cug,  with  the  quantity  a  as  a  parameter.  Eqmtions  (B2) 
should  therefore  be  solved  simultaneously  in  order  to  determine  the 
maximum  obtainable  damping  since 


^1/2 


-0.693  'b 

a  V 


A  convenient  method  for  determining  the  maximum  obtainable  damping  is 
as  follows:  For  several  values  of  a  for  which  the  loops  are  known  to 

be  small,  determine  the  values  of  (Bs  for  which  =  0.  Substitute 

these  combinations  of  a  and  o^g  into  the  expression  ^  and  evaluate 

the  fiinction.  A  plot  may  then  be  made  of  against  a  as  abscissa, 

and  the  value  of  a  for  which  &  =  0  is  the  required  solution.  As 

j  4-  *1^ 

a  check,  this  value  of  a  shoiild  be  substituted  into  —  =  0  and  the 

s 

corresponding  value  of  ojg  determined.  This  combination  of  a  and  cog, 

dTo  - 

if  a  solution,  will  then  satisfy  =  O. 

The  same  analysis  is  applicable  to  the  loops  corresponding  to 
m  =  2,  3  ..  .  .  .  Equations  (Bl)  show  that  R2  and  @2  sre  independent 

of  m:  hence,  ^  is  also  independent  of  m.  The  general  expression 

'  '  (tog 

_  ^  A  /-I1 •4' “t  C  (tkoW 


=  0  is  (see  equations  (B2)) 


2x011  -  Bo 


s  (to. 


Rearranging  the  equation  gives 


2Tan  =  @2 
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dRp 

Thus,  if  the  combinations  of  a  and  0)3  which  satisfy  are  deter- 

s 

mined  for  a  sufficiently  vide  range  of  a  and  substituted  into  the 
right-hand  side  of  equation  (B3)>  the  value  of  a  for  which  loops  will 
cease  to  exist  for  any  value  of  m  can  be  determined.  For  example,  the 


m  =  1  loops  will  break  down  when  02  "  ^s 

when  this  expression  equals  hr:;  and  so  on. 
for  which  the  loops  break  down  decreases  as 


the  m  =  2  loops 
s 

In  general,  the  value  of 
m  increases* 


a 


The  preceding  analysis  was  applied  to  the  m  =  2  and  m  =  3  case 
for  T1/2  =  0.70  second  since  these  loops  were  known  to  have  broken 
down  for  this  damping  (see  fig.  1(e)).  As  was  expected,  the  analysis 
verified  the  fact  that  the  loops  did  not  exist  for  these  values  of  m 
for  T3_^2  =  0.70  second. 
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TABLE  I 

MASS  CHARACTERISTICS  AND  OTABILITY  DERIVATIVES  OF  A 


W/S,  It/sq  ft 
S,  sq  ft  .  . 

t,  ft  .  .  .  . 

p,  slugs/cu  ft 
Y,  ft/sec  .  . 
y,  deg  .  .  . 

Cl  •  •  •  #  • 

•  •  •  •  • 

•  •  •  •  • 

Kxz . 


r\,  deg  . 

•  •  • 

per 

radian 

per 

radian 

per 

radian 

^nr  ^ 

per 

radian 

per 

radian 

per 

radian 

CYp^ 

per 

radian 

Cnp^ 

per 

radian 

per 

radian 

Cng  }  Psr  radian 


TYPICAL  PRESENT-DAY  AIRPLANE 
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Airplane  motion  in  yaw  subsequent  to  an  initial  displacement 
yaw  for  several  combinations  of  k  and  t. 
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Acceleration  autopilot 


Velocity  autopilot 


Angular  frequency,  Ug 


Displacement  autopilot 

Figure  4.-  Types  of  frequency  response  for  |w|  showing  possible  critical 
frequencies  for  acceleration,  velocity,  and  displacement  yaw  autopilots 
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